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On Tightness of Mutual Dependence Upperbound 
for Secret-key Capacity of Multiple Terminals 

Chung Chan 



Abstract — Csiszar and Narayan[3] defined the notion of secret- 
key capacity for multiple terminals, characterized it as a linear 
program with Slepian-Wolf constraints of the related source 
coding problem of communication for omniscience, and upper 
bounded it by some information divergence expression from the 
joint to the product distribution of the private observations. This 
paper proves that the bound is tight for the important case 
when all users are active, using the polymatroidal structure[6] 
underlying the source coding problem. When some users are not 
active, the bound may not be tight. This paper gives a counter- 
example in which 3 out of the 6 terminals are active. 
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I. Introduction 

Consider a group of people with access to correlated dis- 
crete memoryless sources. The problem of communication 
for omniscience[3] (CO) asks how much authenticated pubUc 
discussion is needed so that the active users in the group can 
recover all the sources reliably. The problem of secret-key 
generation[3] (SK) asks how much secret they can agree on. 

Intuitively, the more correlated the sources are, the less 
public discussion is needed for CO because active users 
can learn the sources with their private observations as side 
information. With less information revealed in public, the 
active users can also share more secret. Are the maximum 
savings in public discussion and the amount of secret exactly 
the mutual dependence of the sources? How should one define 
such mutual dependence? 

In [3], Csiszar and Narayan show that the secret-key ca- 
pacity is equal to the maximum savings in public discussion, 
i.e. entropy rate of the sources minus the smallest rate of 
CO. While the active users need not attain omniscience before 
agreeing on a secret key, doing so guarantees certain universal- 
ity result without loss of optimality. They then upper bound 
the capacity by some normalized information divergence [3, 
(26)] from the joint to the product distribution, which is com- 
monly interpreted as mutual dependence of a set of random 
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variables. [6] As will be shown in the sequel, this mutual 
dependence upperbound is indeed tight for the case when all 
users are active. This gives an affirmative answer to the first 
question, and confirms the information divergence expression 
as a measure of mutual dependence. The equivalence is of 
theoretical interest and can simplify the exact computation of 
secret-key capacity in certain special cases as described in [3]. 

In the other case when some users are helpers, it is straight- 
forward to show that the bound is tight for the 3-terminal case 
by an exhaustive test. ' We have found, however, that the bound 
is lose for a 6-terminal counter-example with 3 active users. It 
is unclear if there exists a more general and ubiquitous mutual 
dependence expression that covers this case and other related 
problems. 

II. Preliminaries 

Let Ai := {l,...,m} with m > 2 denote the set 
of all terminals and X^vi := (Xi,...,X„i) be the discrete 
multiple memoryless sources (DMMS), which is a sequence 
of random variables taking values from finite sets. The subset 
A C A4 : 1^1 > 2 denotes the set of active users while its 
complement A'^ is the set of helpers. From [3], the set of 
(strongly) achievable rate tuples of CO for A is, 

^{A) = {Rm e : \/B G ^(^), SW {Rm,13) > 0} 



where 



h{B) 



Each element in ^{A) corresponds to a Slepian-Wolf con- 
straint of CO for A. SW {Rm, 13) is the constraint function 
on the rate tuple Rm {Ri, ■ • • , Rm)- h is the conditional 
entropy function of the DMMS. The Slepian-Wolf constraint 
is satisfied/tight/slack if SW {Rm,B) is >/=/> zero. 

From the set of achievable rate tuples, two key quantities of 
interest, namely the smallest CO rate Rco (A) and secret-key 
capacity Csk{A), can be computed as follows. 



(la) 
(lb) 



Rco{A)= min Vi?, 

l — l 

Csk{A) ^ h{M) - RcoiA) 



'An exhaustive test by a computer program suggests informally that the 
bound is also tight for the 4-terminal case. The 5-terminal case is intractible. 
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Written in matrix form[3], i?co(-4) is the solution to the 
following primal and dual linear programs by the Duality 
Theorem[4, p. 130- 132]. 



(2a) 
(2b) 



Rco{A) 



mm cx 

x:Ax>b 

max yb 

y>0:yA.—c 



(primal) 
(dual) 



where A is an l-hy-m incidence matrix of the Slepian-Wolf 
constraints;- b is the corresponding Z-by-1 column vector of 
conditional entropies; c is a l-by-m row vector of all I's; x 
is a 1-by-TO rate vector satisfying the Slepian-Wolf constraints 
Ax > b; and y is a 1-by-/ non-negative weight vector 

Of particular interest is the following mutual dependence 
upperbound[3, (26)] on Csk(-4). 



(3) I{A) 
where 

/(Ci,.. 



mill HCit- 

2<k<\A\,(Ci,...fik)(^S^k{A) 



,Cfc) := 



k - 
1 



—AY.H{Xc.)-H{Xm) 



,4 = 1 

■D ( P, 



k-i 

^k{A) is the set of all fc-partitionings of Ai such that each 
partition intersects A. I{Ci,...,Ck) is the mutual depen- 
dence with respect to the fc-partitioning (Ci, . . . ,Cfc). P'(-||-) 
is the information divergence from the joint to the product 
distribution. [2] It is the well-known shannon's mutual informa- 
tion in the 2-terminal case, and therefore commonly interpreted 
as mutual dependence for the multi-terminal case. [6] The goal 
of this paper is to confirm this heuristic interpretation by 
studying the tightness of the mutual dependence upperbound. 

III. Statement of results 

From the derivation of the mutual dependence upper- 
bound (3), it is straightforward to see that tightness of the 
bound requires certain Slepian-Wolf constraints to be tight. 
The necessary and sufficient condition is as follows. 

Proposition! (Tightness condition). Csk(^) = I (A) iff there 
exists (Ci, . . . ,Cfc) e 0^k{A) with 2 < fc < |^| and Rm S 
3i{A) such that SW{RM,Cf) = for all i e {1, . . . , k}. 

Not only is this condition useful for a general exhaustive 
test for tightness, it is also fundamental for the subsequent 
proof of the main theorem. The idea is to look for a valid 
partitioning {Ci, . . . ,Ck) such that Csk{A) = I{Ci, . . . ,Ck)- 
This happens iff Cf 's are tight Slepian-Wolf constraints. 

Now, suppose Bi and B2 are two tight Slepian-Wolf con- 
straints. Are their union Bi U B2 and intersection' Bi fl B2 
tight? The answer is affirmative as stated below due to the 
polymatroidal structure[6] . 

is the total number \.^{A)\ of constraints, which is 2™ - 2"'~l-^l - 1. 
^The tightness of the intersection is not essential in the proof of the main 
theorem. 



Proposition 2 (Polymatroidal structure). For all Rj^ e ^(A) 
and BiU B2 €z S§{A), we have 



SW(P^,Si) =0 
SW(P^,S2) =0 



('SW(P^,SiUZ?2) =0 

|sw(p^,z?inz?2) = 



When all users are active, the Duality Theorem and the 
induced tightness of the union of tight constraints implies 
the existence of the desired partitioning for tightness. This, 
however, does not extend to the case when some users are 
helpers. The bound is indeed loose for a particular 6-terminal 
counter-example. Hence, we have the following main theorem. 

Theorem 1 (Tightness). When all users are active, we have 

Csk{M)=I{M) 

When some users are helpers, i.e. A C M., there is a counter- 
example for which Cs]^{A) < I (A) with strict inequality. 

IV. Interpretation of mutual dependence 

Why is a measure of mutual dependence? In the 2- 

terminal case, it is simply the well-known Shannon's mutual 
information /(Xi;X2) := P»(PxiX2 II^^Xi oPx2).[2] But in the 
multi-terminal case, how can we interpret the normalization 
factor and the minimization over different partitionings? 

Consider the sunflower example where X; = (Y, Z^) for i G 
M. The core Y and the petals Zj's are mutually independent. 
For any fc-partitioning (Ci, . . . ,Cfc), 



/(Ci 



,Ck) = 




since H{Xc) = H{Y) + J2,ec ^i^j) for C ^ M,Ci, . . . , Ck- 
Thus, I{A4) = HiY) as expected. From this example, we see 
that the factor l/(fc — 1) compensates for the over-counting 
of H{y) in the sum X)i=i -^(^cj^ and that the optimal 
partitioning need not be unique. 

Suppose (Ci,...,Cfe) is an optimal partitioning that 
achieves the minimum I{A4). Can we say that the random 
variables within the same partition are more correlated than 
those in different partitions? With the tightness result, we can 
recast the question in the CO setting as follows. Consider 
grouping users according to an optimal partitioning. If there 
is a private discussion that leads to omniscience within each 
group before the public discussion, i.e. user j in Ci knows the 
conglomerated sequence of X^, does it reduce the smallest 
CO rate? No because ({!}, . . . , {fc}) is a valid partitioning 
that achieves the same mutual dependence of the set of 
conglomerated random variables. This means that attaining 
omniscience across different groups is the bottleneck of the 
CO problem. One does not need addition redundancy in rate 
to attain omniscience within each group, i.e. SW(PjVi, Cf ) ~ 
for alH £ {1, . . . , fc} and some optimal Rm G ■^{M). This 
gives an operational meaning of the optimal partitioning. 

V. Counter-example 

In this section, we describe a counter-example for which 
the mutual dependence upper bound Csk{A) < I (A) is 
loose, i.e. satisfied with strict inequality, when some users are 
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helpers. Let m = 6, ^ = {1,2,3} and (Yi, Y2, Y3, Y4) be iid 
uniformly random bits. The DMMS is defined below as 
the XOR © of every distinct pair of the random bits. 



Proof: Consider a partitioning (Ci, . . . ,Cfc) £ ^fe(^). 



Xi 


- Yi 




©Y3 




X2 


= Yi 






©Y4 


X3 






Ya 


©Y4 


X4 




Y2 


©Y3 




X5 




Y2 




©Y4 


Xe 


= Yi 


©Y2 







There are altogether / ~ 55 Slepian-Wolf constraints B G 
S§{A). With some algebra, the conditional entropy function 
h can be simplified to the following expression that depends 
only on the cardinality of the constraint. 



(4) 



if|S|e{l,2} 
h{B) = {l if|S|e{3,4} 
2 if \B\ = 5 



Solving the linear program (la) and applying (lb), 

9 _ 
4 ~ 

which is achieved by the unique optimal rate tuple with 



Rco{A) 



Csk{A) = 1.75 



Ri — R2 — — 



R4 = R5 — Rq — - 



Appendix B describes how to compute and verify this solution. 

It remains to show that I{A) > 1.75. Consider a 3- 
partitioning (Ci,C2,C3) G ^^^{A) such that each partition has 
2 elements. Applying (4) with \Cf\ = 4, 



IiCi,C2,C3) = h{M)- 

. 1 + 1 



1 



3- 
2.5 



1 



All other cases can be computed analogously as follows by 
enumerating the integer partitionings of 6. 



integer partitioning {|CJ|, . . . , |C^|} I{Ci 



2 + 2 + 2 


{4,4,4} 


4- 


1+1+1 _ 
2 


= 2.5 


3 + 2 + 1 


{3,4,5} 


4- 


1+1+2 
2 


= 2 


4+1 + 1 


{2,5,5} 


4- 


0+2+2 


= 2 


3 + 3 


{3,3} 


4- 


(1 + 1) 


= 2 


4 + 2 


{2,4} 


4- 


(0+1) 


= 3 


5 + 1 


{1,5} 


4- 


(0 + 2) 


= 2 



We have the desired strict inequality that, 

liA) = 2 > 1.75 = Csk{A) 

Appendix B gives an alternative argument that explains how 
this counter-example is constructed. 

VI. Proofs of results 

The mutual dependence upperbound is derived in [3] by 
removing certain Slepian-Wolf constraints as follows. 

Proposition 3 (Mutual dependence upperbound). 

(5) Csk{A) < I{A) VAciM:\A\>2 



i=i jecf 



E E 

i=i iJeM jec, 
(fc - 1) ^ R, 



(6) 



3&M 



Applying only the Slepian-Wolf constraints Cf s on Rm to 
the R.H.S. gives the desired bound. ■ 
The tightness condition then follows from the tightness of 
those Slepian-Wolf constraints used to derive the bound. 

Proof of Proposition 1: Consider proving the 'if case. 
Let (Ci,...,Cfe) e ^k{A) and Rm & ^{A) be such that 
SW(i?^, Cf) = for alH e {1, . . . , k]. Then, from (6), 



j=i 



Rj=h{M)-I{Ci,...,Ck) 



< h{M)-IiA) 

By definition RcoiA) < J2Zi < HM)-I{A). Together 
with (5), we have Rco{A) = H{Xm) - HA). 

Consider proving the contrapositive of the 'only if case. 
Suppose for all 2 < fc < |^|, (Ci, . . . ,Cfc) G ^k{A), Rm e 
^{A), there exists i e{l,...,k} such that SW{Rm , Cf) > 0. 
Then, for all optimal Rm, 



RcoiA) 



j=i 



Rj > H{Xm)-I{A) 



which implies Rco{A) ^ H{Xm) ~ I (A) as desired. ■ 
The polymatroidal structure[6] of the source coding problem 
can be used to prove that union and intersection of tight 
Slepian-Wolf constraints are tight as follows. 

Proof of Proposition 2: Consider some Rm G ^{A) 
and BiU B2 G MiA) such that SW{Rm,Bi) = and 
SW{Rm,132) = 0. Then, Bi U B2 G ^{A) implies that 
Bi U B2, Bi, B2 and Bi D B2 do not contains A. If they 
are non-empty, they are all in ^^{A) by definition. Let 
r{B) := E,;es^» for all B C M. Then, 

r{Bi)^h{Bi) r{BiUB2)>hiBiUB2) 
riB2) = KB2) r{Bi n B2) > h{Bi H B2) 

To show that the union is a tight constraint, 

h{Bi U B2) < r{Bi U B2) - r{Bi) + r{B2) - r{Bi n B2) 

<h{Bi) + h{B2)-h{Bir\B2) 
< h{Bi U B2) 

where the last inequality is by the supermodularity of h. (see 
Lemma 1 in Appendix A) Hence, r{Bi U S2) = h{Bi U S2). 
Similarly, r (Si 0^2) = h{BlC^B2)■ The case when Si 0^2 = 
holds trivially. ■ 
The idea of the proof for the main theorem is to first obtain 
an initial set of tight Slepian-Wolf constraints from the dual 
linear program. Then, use the fact that union of tight Slepian- 
Wolf constraints are tight to construct the desired partitioning 
for the tightness condition. 
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Proof of Theorem 1: Consider some optimal solution y 
to the dual linear program (2b) with A = Ai. Let < > 1 be 
the number of non-zero entries in y. Construct the 1-by-^ row 
vector y by eliminating the zero entries in y; and the i-by-m 
submatrix A of A by eliminating the corresponding rows, i.e. 
removing the i-th row if = for i e {1, . . . , ^}. 

Let A,j be the j-th column of A. Construct the desired 
partitioning (Ci,...,Cfc) £ ^k{M) by partitioning A into 
classes of identical columns and let the i-by-1 vector sj be 
the column of the i-th class. The precise construction is as 
follows. 

Input: A 

Output: k, (si, . . . , Sfc) and (Ci, ...,Ck) 
k := 1; Ci := 0; Si := l^^; 
for j := 1 to m do 

if A,j = sJ for some i < k then 
I Add j to C; 
else 

I fc := /c + 1; Cfc := 0; := A^f, 
end 

end 

To argue that fc > 2, note that A does not have any rows of 
all I's nor rows of all O's because A4 and are not Slepian- 
Wolf constraints. With i > 1, at least two columns of A are 
distinct. The other constraint that k < m holds trivially. 

It remains to argue that Cf s are tight Slepian-Wolf con- 
straints in essence of the tightness condition in Proposition 1 . 
As an immediate consequence of the Duality Theorem (see 
Lemma 2 in Appendix A), rows of A correspond to tight 
Slepian-Wolf constraints. Since unions of tight constraints are 
tight by the polymatroidal structure (see Proposition 2), it suf- 
fices to show that Cf 's are unions of constraints corresponding 
to rows of A. 

To do so, define ^ as elementwise negation and o as the 
logical matrix multiplication, in which addition and multipli- 
cation are replaced by logical 'or' and 'and'. Then, (^Si)o A 
corresponds to taking union of constraints in A whose corre- 
sponding entry in Si is 0. It suffices to show that s^o A indeed 
corresponds to Cf, or equivalently, that (-is^/)© sJ — 1. 

Assume to the contrary that there exists i i' in {1, . . . , k} 
such that (-iSi')© sJ = instead. Since sj's are constructed 
from distinct columns of A, there exists j ^ j' in Ai such 
that, 

(-iA,j')'^s A,j = but A,j> ^ A,j 

This implies that the constraint of A,j' is a proper superset 
of that of A,j. Since y > 0, 

yA,j, > yA.j 

But this contradicts yA = yA = 1, which is the constraint 
of the dual linear program. This completes the proof. (See 
Section V for the counter-example.) ■ 

VII. Conclusion 

The mutual dependence upperbound I{A4) on the secret- 
key capacity in [3] is proved to be tight for the case when 



all users are active. This gives an operational meaning to 
the mutual dependence expression, and therefore confirms its 
heuristic interpretation as a measure of correlations among a 
set of random variables. 

The proposed proof uses the polymatroidal structure in the 
source coding problem pointed out by [6]. Starting with an 
arbitrary solution to the dual problem (by the Duality Theorem 
in linear programming already mentioned in [3]), an initial 
set of tight Slepian-Wolf constraints is obtained. The desired 
set of tight Slepian-Wolf constraints is then derived using the 
polymatroidal structure. 

As shown by the counter-example with 3 active users and 3 
helpers, the mutual dependence upperbound need not be tight 
for the case when some users are helpers. Thus, the mutual 
dependence expression I{A) in this case is not supported 
with the operational meaning of CO and SK, even though 
an exhaustive test shows that it is tight for the 3-terminal case 
and (informally with the help of a computer) for the 4-terminal 
case. It is unclear if there is a more general mutual dependence 
expression that cover this case or other problems related to the 
mutual dependence of a set of random variables. 

Appendix A 
Clarifications 

Lemma 1 (Supermodularity). For all 81,82 C A4, 
h{8i) + h{82) < h{8i U 82) + h{8i n 82) 
Proof: Subtracting L.H.S. from R.H.S. gives, 

h{8i U 82) + h{8i n 82) - h{8i) - h{82) 
= -i/(Xsjn6^) - iflXgjus^) + H{Xb^) + HiXei) 

-i/(XBjuB^|Xe5;nB^) + -H^(Xi3j IXgjnei) + H {XBi\^BinBi) 
— ^(Xbj ; Xgc jX^jne^ ) 

which is positive as desired. ■ 

Lemma 2 (Complementary slackness[4, p. 135- 136]). Con- 
sider the primal (2a) and dual {2h) formulations for Rco{A). 
For all optimal solutions x and y, i € {1, . . . , if yi > 0, 
then the i-th row Ai, of A corresponds to a tight Slepian-Wolf 
constraint, i.e. Ai,x — bi. 

Proof: By the (strong) Duality Theorem, cx = yb. Since 
yA = c for y to be feasible, we have yAx = yb, or 
equivalently, 

I I 
^yjAj,x ^^Vjbj 

Assume to the contrary that there exists i E {!,...,/} such 
that yi > and Ai,x ^ hi. Since Ax > b, and y > 0, 
the L.H.S. of the last equation would instead be strictly larger 
than the R.H.S., which is a contradiction. ■ 

Appendix B 
Computations for the counter-example 

In this section, we compute i?co(-4) for the counter- 
example in Section V using the Multi-Parametric Toolbox[7] 
for Matlab. First, we initialize the toolbox to use an LP solver 
that always terminates at a vertex optimal solution. 
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% initialize tiie LP solver with CDD Criss-Cross 
opt=mpt_lnit ( ' ipsoiver' , ' cdd' ,'abs_tol' ,le— 8); 

The matrices for the linear program (2a) can be constructed 
manually as follows. 



I 














();() 


I 











0;1 


1 




















1 








0;1 





1 














1 


1 








0;0 








1 





0;1 








1 











1 





1 





0;1 


1 





1 














1 


1 





0; 1 





1 


1 





0;0 


1 


1 


1 




















1 


0;1 











1 








I 








1 


0;1 


1 








1 


0;0 





1 





1 





1 





1 





1 


();0 


1 


1 





1 














1 


1 


0;1 








1 


1 


0;0 


1 





1 


1 





1 


1 





1 


1 


0;0 





1 


1 


1 





1 





1 


1 


1 


0;0 


1 


1 


1 


1 


0;0 














1 


1 














i;0 


1 













I 


I 











1;0 





1 








i;l 





1 








1 





1 


1 








1;0 








1 







I 








1 





l;0 


1 





1 





1;1 


1 





1 





1 








1 


1 





i;i 





1 


1 










1 


1 


1 





1;0 











1 


i;i 











1 


1 





1 








1 


i;i 


1 








1 










1 





1 


1;1 





1 





1 


1;0 


1 


1 





1 


1 











1 


1 


i;i 








1 


1 







1 





1 


1 


1 ; 1 


1 





1 


1 


i;0 





1 


1 


1 


1 


1 





1 


1 


1 


i;0 


1 


1 


1 


1 





00000000000 1 0000000 1 000 1 0,... 

1 1 I 1 1 1 1 1 1 1 1 1 2 1 2 2 



c = [l 11111]; 

We can then use the command mpt_solveLP to solve the 
primal and dual linear programs. To match the API, y is 
a column vector instead of a row vector, and nRCO is the 
negation of Rco{-^)- 

% Primal: Rco — miria; c * x subject to: —A * x <— —b 
[x , Rco] = mpt_solveLP (c ', — A,— b ) ; 

% Dual: nRco — miriy b' ^ y subject to: A' ^ y — c' 
lb = zeros(size(A,l) ,1); 
ub=lnf*ones(size(A,l) ,1); 

[y ,nRco] = mpt_solveLP(-b ,[],[], A' , c ',[],[], lb , ub ) ; 

The optimal solution x== [.25 .25 .25 .5 .5 .5] of 
the primal achieves Rco==2.2 5. The optimal solution of 
the dual has non-zero entries equal to 0.25 at positions 
of the tight Slepian-Wolf constraints. One can check that 
y (A*x-b<=opt . abs_tol) returns a sequence of all 0.25's, 
and sum (A*x-b>opt . abs_tol) is 0. Now, to verify that 
the computed x is indeed optimal (since the output from a 
computer cannot be used as a formal proof), one can check that 
X and y are both feasible and satisfy the primal/dual optimality 
criteria in Theorem 2.9 of [5, p.48]. This implies that both x 
and y are optimal, and i?co(-4) = 2.25 as desired. The mutual 
dependence upperbound computed in Section V is therefore 
formally proved to be loose. 

There is an alternative explanation that the bound is loose 
without calculating it explicitly. The underlying reasoning has 
guided the construction of the counter-example, and the proof 
of tightness in the case when all users are active. Consider 
the tightness condition in Proposition 1. If one can show that 
for all optimal solutions in 3i{A), there is no subset of 2 < 
k < \A\ tight Slepian-Wolf constraints whose complement 
partitions A4, then the mutual dependence bound is loose. 
For the particular optimal solution x, exactly six Slepian-Wolf 
constraints are tight, i.e. 







A:-- 








x: — 




b: = 


"1 





1 


1 





0" 




-1- 

4 




"l1 





1 


1 





1 







1 

! 
! 




1 


1 
1 


1 
1 








1 




1 


1 
1 






1 

2 


1 





1 


1 


1 


1 




? 
? 

.2. 




2 





1 


1 


1 


1 


1 
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where x is the optimal solution x; Rows of A corresponds 
to the tight Slepian-Wolf constraints; and elements of b are 
the corresponding conditional entropies. Note that rows of A, 
has either one or three O's. For the complements of a subset 
of 2 or 3 Slepian-Wolf constraints to partition M, the only 
possibility is to have two tight Slepian-Wolf constraints whose 
incidence vectors have three O's. But there is no two such rows 
of A having O's at complementing positions. Thus, the bound 
is loose if x is the unique optimal solution, which can be 
shown using the PUFAS algorithm in [1] as follows. 

We first express the primal linear program in equational 
form by introducing a slack variable for each constraint to 
take up the slack. Let Xs be the 55-by-l column vector of 
slack variables. Then, 



RcoiA) 



mm cx 

Ax>b 



min [c 







where I is an identity matrix with matching dimensions. The 
vertex optimal solution xl of this equational form can be 
computed as follows. 

% Solving the primal problem in equational form 
ceq = [c zeros ( 1 , size (A, 1 ))] ; 
Aeq=[A — eye(size(A,l))]; 
lbeq=zeros(size (Acq ,2) , 1 ) ; 
ubeq=Inf*ones(size (Acq , 2 ) , 1 ) ; 

[ xl , Rcol ]=mpt_solveLP(ceq ' , [ ] , [ ] , Acq , b ,[],[], Ibeq , ubeq ) ; 

One can verify that Rcol==2.25, xl(l:6)==x and 
xl (7 : 61) ==A*x-b within the absolute tolerance. If the so- 
lution is not unique, there will be a different vertex optimal 
solution with at least one positive entry at the position where 
the corresponding element of xl is zero. Let d be the incidence 
row vector of O's in xl. Maximizing dx subject to an 
additional constraint for optimality that [ c o ] [ equals Rcol 
gives an alternative solution if there is one. If the solution is 
unique, the optimal value would be 0, achieved by xl. This 
test can be implemented as follows. 

% Verifying uniqueness of optimal solution x (i.e. show notUnique=0) 
[x2 , no tU nique ] = mpt_solveLP(— double ( xKle — 8 ),[],[],.. . 
[ Aeq ; ceq ] , [ b ; Rcol ],[],[], Ibeq, ubeq ) ; 

One can check that notunique==0 and xl==x2 within 
the absolute tolerance. This shows x is the unique optimal 
solution, and hence the mutual dependence bound is loose. 
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